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ABSTRACT 



in 

O 

. After Michel (1994) introduced a phenomenological picture of "rapid 

^ ! magnetization" of newly born neutron stars (NSs), Muslimov & Page (1995) 

^! suggested that the physical conditions accompanying the formation of a NS 

2 ! (convection in the NS envelope, high-rate accretion in a supernova) may 

^ ! result in the large scale surface magnetic field of the NS having a low value 

^ ' (~ 10* — 10^ G), while the bulk of the magnetic flux is submerged under the 

stellar surface. The subsequent growth of the large scale surface magnetic field 
^ ! then occurs owing to ohmic diffusion of a strong internal field (~ 10^^ — 10^^ G) 

I and proceeds over a timescale of 10^ — 10^ yr, depending on the early thermal 

history of the NS, initial distribution of magnetic flux, and electrical resistivity 
in the crust. Referring to this suggestion, we perform numerical computations 
to demonstrate its possible relevance to young radio pulsars. In our calculations, 
we use different NS structures based on a model of dense matter presenting 
a phase transition to kaon condensation which softens the equation of state. 
This model implies both slow and fast neutrino cooling, depending on the 
stellar mass (which is assumed to be in the range of 1.35 — 1.50 Mq). We 
present sequences of secular evolution of the surface magnetic field, spin-down 
luminosity, spin period, "braking index," and spin-down age best matching the 
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corresponding measured and derived quantities for PSR 0531+21 (Crab), PSR 
0540-69, and PSR 1509-58. We conclude that the effect under consideration 
reproduces remarkably well both the observed rotational characteristics and 
derived values of the surface magnetic field strength for these radio pulsars, the 
only ones having reliable measurements of the braking indices. In addition, our 
analysis suggests that NSs in the Crab and PSR 0540-69 have experienced fast 
neutrino coohng and that their masses are above the critical mass for the phase 
transition while the NS in PSR 1509-58 has a mass below this critical mass and 
has therefore undergone slow neutrino cooling. 

Subject headings: pulsars: individual (PSR 0531+21, PSR 0540-69, PSR 
1509-58) — stars: neutron — stars: magnetic fields 
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1. INTRODUCTION 

The timing of radio pulsars reveals that their "clock mechanisms" are subject to 
a systematic delay (spinning down) and that various irregularities in their run such as 
"glitches," "microglitches," etc., may occur. The measurements of the first and second time 
derivatives of pulsars' spin period, P and P, respectively, can provide us with invaluable 
information about the rotational dynamics of neutron stars (NSs). For the standard formula 
describing the secular decrease in the angular velocity Vt of an NS, 

i7oc-l]" , (1) 

where n = 3 corresponds to the mechanism of pure magnetic dipole braking, the "braking 
index" (see Manchester & Taylor 1977, hereafter MT, eq. [6-17]), 

can, in principle, be determined from observations. In practice, the measurements of 
braking indices are contaminated by the effects of "restless" behavior of P over short 
timescales and are a point of controversy. Only three young pulsars have measured values 
of n not dominated by timing noise: PSR 0531+21 (Crab), n = 2.509 ± 0.001 (Lyne, 
Pritchard, & Smith 1988); PSR 0540-69, n = 2.01 ± 0.02 (Nagase et al. 1990; Gouiffes, 
Finley, & Ogelman 1992); and PSR 1509-58, n = 2.837 ± 0.001 (Kaspi et al. 1994). The 
observed rotational properties of these three pulsars are listed in Table |I]. 

Theoretically, as is summarized in MT (see Table 9-1), there are a number of factors 
affecting the braking index, e.g., multipole electromagnetic radiation {n > 5), gravitational 
quadrupole radiation (n = 5), magnetic field decay (n > 3), radial deformation of field 
lines [1 < n < 3), pulsar wind {n < 3), relaxation of the NS equilibrium form [n < 3), 
pulsar transverse velocity (n < 3), etc. To these, one should add two mechanisms recently 
proposed: spin-up of the NS during its early evolution owing to the intense neutrino 
emission {n < 0)(Alpar & Ogelman 1990) and crustal "plate" motions (Ruderman 1991), 
the latter being able to produce values of n larger or smaller than 3. Also, Blandford & 
Romani (1988) have discussed the possibility that the braking indices of PSR 0531+21 
and PSR 1509-58 may be consistent with a secular growth of their magnetic field. They 
appealed to the thermomagnetic instability in the NS crust (Blandford, Applegate, & 
Hernquist 1983, see also Urpin, Levshakov, & Yakovlev 1986) as a basic reason for the 
magnetic field growth. Finally, Michel (1991) has also emphasized that the monotonic 
increase of the magnetic field with time would give braking indices less than 3 in all young 
pulsars where it can be determined. 
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In this paper we employ our recent calculations of growth of the surface magnetic 
field in a newly born NS (Muslimov & Page 1995, hereafter MP) to demonstrate the 
striking agreement of our model with the available observational data for the Crab, PSR 
0540-69, and PSR 1509-58. Our tentative analysis (see MP) had been stimulated by a 
phenomenological picture recently discussed by Michel (1994, hereafter M94) that the 
magnetic field of a NS at birth is low and increases, by a factor of ~ 10^, to the typical 
value of ~ 10^^ G during a period of a few hundred years following the supernova. MP 
suggested, as a physical basis for this picture, that the surface magnetic field of a young 
NS may grow from a very low value (~ 10^ — 10^ G) owing to ohmic diffusion of a strong 
internal magnetic field (~ 10^^ G) that was initially submerged under the surface layers. 
Assuming that an NS spins down as a result of magnetic dipole radiation, we determine the 
free fitting parameters of our model. For example, for each of these pulsars we specify the 
initial spin period and initial profile of the magnetic flux in the crust. 

These early magnetic and spin evolutions are strongly affected by the thermal evolution 
of the NS. In the present study we exploit different stellar structures and thermal histories 
to fit the data for the three young pulsars. Specifically, we consider different NS masses and 
the kaon condensation phase transition (Thorsson, Prakash & Lattimer 1994; hereafter 
TPL) in the NS core which results in different cooling histories depending on the mass. 
Within the framework of our analysis, the fast cooling of an NS provides the best fit for the 
Crab and PSR 0540-69, including accounting for the discrepancy between their "real" and 
spin-down ages. In the case of PSR 1509-58, we will need a slow coohng, i.e., a star of mass 
below the critical mass for the phase transition. 

In § H we present our numerical model, and in § |^ we illustrate matching of our 
calculations with the observed parameters for PSR 0531+21 (Crab), PSR 0540-69, and PSR 
1509-58. Our principal conclusions are summarized in § ^. 

2. BASIC ASSUMPTIONS AND METHOD 

2.1. General Picture 

We assume that an NS is born with a spin period in the range 20-40 ms. Alternatively 
(see, e.g., M94), the luminosities of the NSs born as rapid rotators will become large 
before the NSs slow down to the present luminosities of the Crab, PSR 0540-69, and PSR 
1509-58 which are already large enough. To calculate the spin evolution of the NS we use 
the standard formula (see, e.g., MT, p. 110, 180) for the spin-down torque resulting from 
emission of magnetic dipole radiation. 



If the secular evolution of the dipole magnetic moment of an NS is the only physical 
reason affecting its spinning down, it is straightforward to derive the following expression 
for the braking index 



where Bgurf is the surface value of the magnetic field strength at the magnetic pole. (A 
similar formula, but missing a factor 2 in the second term, has been used by Chanmugam 
& Sang [1989] in their analysis of the possible response of braking indices to ohmic field 
decay in pulsars.) 

In formula (0) Cl is always negative and, depending on whether Bsurf > (field growth) 
or Bsurf < (field decay), we have n < 3 or > 3, respectively. Here we attempt to 
quantitatively examine the possibility that, in a newly born NS, the magnetic field (of 
strength ~ 10^^ — 10^'^ G), initially trapped under the stellar surface, diffuses ohmically up 
through the crust ( MP). The main uncertainty in this picture is the form of the initial 
distribution of the magnetic fiux in the upper crust. For illustrative purposes, we describe 
the initial profile of magnetic fiux in the crust by referring to the depth down to which the 
fiux is submerged under the stellar surface. Approximately at this depth the magnetic fiux 
has a maximum value and decreases by a factor of ~ 10^ — 10^ toward the stellar surface. 



We now consider the evolution of the dipole component of a purely poloidal magnetic 
field B in the NS crust. For a dipole mode, in spherical coordinates (r, 6, 0), we can express 
B in the form 



Zr or J 

where Bq is some normalization magnitude of the field, and are unit vectors in the 
radial and meridional directions, and S = S{r, t) is the stream function. 

Using the standard equation for the magnetic field evolution in fiat spacetime 
(neglecting the hydrodynamic motions, thermomagnetic effects, and anisotropy of the 
electrical conductivity) , 




(3) 



2.2. Basic Equations and Solution 




(4) 




(5) 



we arrive at the following equation for the evolution of S: 
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where a = o"(r, t) is the electrical conductivity of matter. 

We solve equation (|^) numerically subject to the appropriate boundary conditions at 
the stellar surface and at the base of the crust. At the surface (r = R) we impose the 
standard boundary condition that the internal field merges continuously with an external 
vacuum field. For a dipole field this condition reads 

R +S = . (7) 

or 

At the base of the crust (r = Vb) we assume that 

5 = Const . (8) 

This second condition is perfectly justified for our problem, and it means that the bulk of 
the magnetic flux [the magnetic flux through a hemisphere of radius r is proportional to 
the value of the function S{r,t)] is frozen into the inner crust (outer core), at least over a 
timescale of ~ 10® yr. The ohmic diffusion and decay of the magnetic field in the innermost 
layers of the crust occur on a much longer timescale and are of no importance for our 
present purpose. Note that in our calculations the conductivity of the crust depends on 
both the density and temperature of matter (see, e.g., MP and references therein for a 
relevant discussion of the conductivity regimes in the crust and corresponding numerical 
estimates), the latter changing with time as the NS cools down. 

We perform our calculations of the spin-down luminosity, spin period, braking index, 
and spin-down age (r^^ = P/2P) of the NS by solving self-consistently the system of 
equations describing both its magnetic and spin evolution (owing to magnetic dipole 
braking). We also employ formula (^ in our field evolution code to calculate the secular 
evolution of the braking indices. 

The principal parameters needed to be specified in our computations are the initial 
spin period of the NS, Pq, and the initial surface value of the magnetic field strength, B^^^j-. 
We assume that for the Crab, PSR 0540-69, and PSR 1509-58 the value of 5°^^^ is in the 
range of (1 — 5) x 10^° G, while Pq ranges from ~ 20 to ~ 40 ms. It is important that our 
choice of Pq for the Crab and PSR 0540-69 is not arbitrary, and we shall discuss this below. 

Our calculations show clearly that the character of early magnetic and spin evolutions 
of a NS is determined by the dynamics of freezing and the structure of its crust. For 
example, the stellar models allowing for a rapid cooling result in a power-law temporal 
growth of the surface value of the magnetic field strength in the saturation regime. In 
contrast, the field evolution in the saturation regime for a relatively slowly cooling NS is 
more complex and cannot be described by a simple power-law time dependence. Also, 
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the fact that the NSs in the Crab and probably in PSR 0540-69 have spin-down ages well 
exceeding their real ages can be understood naturally within the framework of our analysis 
only if we assume that these NSs have undergone a stage of rapid neutrino cooling. The 
corresponding stellar models and thermal histories we have employed in our calculations for 
these pulsars do reproduce remarkably well the discrepancy between their spin-down ages 
Tsd and expected real ages tmodei (see § H). Finally, in the case of PSR 1509-58 our results 
may testify that the spin-down age of the associated NS is less than its real age and that 
this can be attributed to relatively slow neutrino cooling of the NS. 

The standard equation describing the evolution of the spin period owing to magnetic 
dipole losses can be easily integrated to give 

P = Po{l + al^Bl,fdt'Y' , (9) 

where a is a numerical factor depending on the radius, moment of inertia, and initial value 
of the spin period of the NS. Using the above expression, we can represent the spin-down 
age of an NS in the form 

P K 



2P iL^io^-^*'- 

where K = [\ — {Pq/ PY]~^ . Note that in our problem the integral Jq B'^^j.^dt' is dominated 
by the contribution from the saturation regime, and the integration can therefore be 
performed over the interval of time during which the field has been saturating. As we have 
emphasized above, for the stellar models with a rapid neutrino cooling (that have used for 
NSs in the Crab and PSR 0540-69) the surface value of the magnetic field strength increases 
in the saturation regime according to a power low. Thus, assuming that in the saturation 
regime, for NSs in the Crab and PSR 0540-69, Bg^rf oc t^, we obtain 

K 

Tsd 1 (11) 

and 2 

n ~ 3 —K . (12) 

2e + 1 ^ ' 



From equations (0) and (|T2D we find that 



/s:~3-n + — (13) 

t 

and 

^~ 2(n-3 + A-) ■ 
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We can use the expressions (|T3|) and (|T^ for rough estimates of the values of Pq and 
e for the Crab and PSR 0540-69. Substituting the values P, n, Tsd, and t = tmodei into 
equations (13) and (14), we find that for these pulsars Pq ~ 23 ms and e ~ 0.18, and 
Pq ~ 41 ms and e ~ 0.27, respectively. These estimates justify our choice of the range of 
the initial spin periods for the Crab and PSR 0540-69. However, the above reasoning is not 
applicable to the NS in PSR 1509-58 because of the needed relatively slow stellar cooling 
and, as a consequence of this, more complex temporal evolution of the surface value of its 
field strength in the saturation regime. In this sense, our choice of Pq for PSR 1509-58 is 
quite arbitrary. 

For a given range of variation of Pq and B^^^p we survey different profiles of the initial 
distribution of the magnetic flux in the crust to fit the observed and derived parameters 
such as P, n, Tsd, and Bgurf- As a result of such a fitting we are able to determine the 
sensible range of the initial depths of submergence {zsub) of the magnetic field in the crust 
(see Table |l]). Note that the rate of spinning down and the value of the braking index 
decrease with increasing the initial depth of submergence Zguh of the flux (for a given NS 
structure and thermal evolution), with the values of braking indices being most sensitive to 
the variation of this depth as a parameter. In our calculations, the initial profile and depth 
of submergence of the magnetic flux in the crust are therefore determined by matching the 
theoretical and measured values of braking indices for a given (derived) spin-down age of 
the pulsar. Note also that the variation of Pq only slightly affects the calculated values of 
the spin period and braking index for relatively long-period pulsars such as PSR 1509-58. 

2.3. Neutron Star Models and Cooling 

Our many trials to fit the observed parameters of the three young pulsars led us to 
focus on different NS thermal histories and stellar structures. To perform a consistent 
analysis, we are thus forced to assume that the three NSs have different masses and that 
the equation of state (EOS) allows fast or slow neutrino cooling depending on the mass. 

A number of emission mechanisms have been proposed for the fast neutrino cooling 
(see, e.g., Pethick 1992, Page 1994). In this paper, we adopt the possibility of kaon 
condensation which is presently the scenario most seriously discussed in the literature 
(see Brown & Rho 1995 for a review). To build a consistent EOS, we follow the method 
of TPL, in which the effect of kaon condensation is superposed on an underlying purely 
baryonic EOS in the mean field approximation (very similar results have been obtained 
with a different formalism by the Kyoto group [see, e.g., Tatsumi 1995]). Due to the large 
uncertainties in the structure of matter at the densities characteristic of NS cores, we have 
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many free choices. We select the baryonic EOS in a standard way, while we parameterize 
the kaon condensation properties to satisfy our needs, namely: 

1. We select, as TPL did, the baryonic EOS from the set of parametric models of 
Prakash, Ainsworth & Lattimer (1988, hereafter PAL): we use the model "PAL33" with 
an incompressibility Kq = 240 MeV and a symmetry energy function F{u) = u^^^ (where 
u = n/uQ, n is the baryon number density and uq = 0.16 fm~^ is the "saturation" density). 
This EOS is intermediate in stiffness between the standard microscopic non relativistic 
FP EOS (Friedman & Pandharipande 1981) and the microscopic relativistic MPA EOS 
(Miither, Prakash & Ainsworth 1987). It provides a proton fraction that is low enough to 
prevent the fast neutrino emission by the direct Urea process (Lattimer et al. 1991). The 
mass-radius relationships for these three EOS are shown in Figure 

2. The properties of the kaon condensate in the formalism developed by TPL depend 
strongly on the values of the coupling constants of the model Lagrangian, all of them 
but one being fortunately constrained by experiments. The only poorly constrained 
parameter is thus (expressed as a^rris, where rris is the mass of the strange quark), and 
its reasonable range of variation is from —134 MeV down to about —300 MeV (see TPL). 
Given the underlying baryonic EOS, the value of asnig determines the critical density for the 
condensation phase transition, i.e., the critical mass for an NS to contain a kaon condensate 
core. Note that the maximum NS mass also depends on a^rus but not so sensitively as the 
critical mass. We adopt in our calculations the value of asuis = —190 MeV. The resulting 
mass-radius relationship and the density profiles for five models of different masses are 
plotted in Figures |l| and ^ The NSs with masses below 1.37 Mq , the critical mass for 
appearance of the condensate, are "standard" NSs, while more massive ones should rather 
be called "nuclear" stars (Brown & Bethe 1994), since they contain roughly equal amounts 
of protons and neutrons, as laboratory nuclei do, and a very strong kaon condensate (unlike 
laboratory nuclei). 

The motivation for our particular choice of a^rris = —190 MeV is twofold: first, it 
provides a transition from the slow to the fast neutrino cooling at a mass of 1.37 Mq , 
which is within the range of measured NS masses (van Kerkwijk, van Paradijs & Zuiderwijk 
1995). Second, it gives a maximum mass of 1.575 Mq that is still consistent with the data 
and is also compatible with the possibility that the NS formed in SN 1987A went into a 
black hole (Bethe & Brown 1995). 

With regard to the cooling we come up with the following scenario. A NS of mass 
M < 1.37 Mq allows only the slow modified Urea process along with the associated nucleon 
bremsstrahlung processes and therefore follows the "standard" cooling scenario . If M > 
1.37 Mq , then the presence of the kaon condensate increases the neutrino emission by 
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orders of magnitude (Brown et al. 1988) resulting in fast cooling (Page & Baron 1990). 
Above 1.463 Mq the direct Urea process is also allowed and dominates the total neutrino 
emissivity that results in an even faster cooling (Page & Applegate 1992). The cooling 
curves for five models are shown in Figure ^ where the occurrence of the three different 
scenarios is clearly seen. The evolution of the profiles of internal temperature for two of 
these models is shown in Figure The profiles have strong gradients during the early 
phases of cooling, but after a few years the crust temperature in the density range 10^^ 
- 10^^ g-cm~'^, where the currents producing the magnetic field are located, is uniform to 
within 50% or better. 

Some details of the method used to produce these cooling models can be found in 
Page & Baron (1990) and Page & Applegate (1992). In this paper we introduce some 
modifications, which, besides the changes in the EOS, are as follows. We use a proper 
treatment of the suppression of the specific heat and neutrino emission by baryon pairing 
according to Levenfish & Yakovlev (1994a,b). For the kaon condensate neutrino emissivity 
we use the formula 

= 7.6 X 10^^ X Tl f (^) ^^[1 + ^{g'^f] sm' Oc erg cm'^ s'\ (15) 

which is obtained from the similar formula of Brown et al. (1988) by replacing the 9\ 
by sin^ 9^ to take into account that 9k is not small. Also, an additional factor of 2 has 
been included following Thorsson et al. (1995). In the above formula m„ and m* are, 
respectively, the neutron free and effective masses; /ig is the electron chemical potential; m^r 
is the pion mass; ~ 1 is the in-medium Gamow- Teller coupling constant; = 13° is the 
Cabbibo angle; Tg is the temperature in 10^ K; and 9k is the chiral angle characterizing 
the kaon condensate strength {9k grows rapidly to near 90° above the critical density). We 
have incorporated nucleon pairing by employing the results of, correspondingly, Ainsworth, 
Wambach, & Pines (1989) [for neutron ^Sq pairing (inner crust)], Takatsuka (1972) [for 
neutron ^P2 pairing (core)], and Takatuska (1973) [for proton ^Sq pairing (core)]. These 
core neutron and proton pairings (see Page 1994 for a comparative presentation of the 
relevant calculations) do not affect our computations because at the densities typical for 
the inner core (where the intense neutrino emission occurs) the corresponding critical 
temperatures Tc are vanishingly low and do not allow the pairings. In contrast, neutron 
pairing in the inner crust is more important, since it reduces substantially the specific heat 
in this region, where the main currents producing the stellar magnetic field are located. As 
a result, the cooling of the crust proceeds faster during the first few decades. Note that the 
calculation of Ainsworth et al. (1989) is the most accurate treatment available to date and 
gives Tc of the order of lO^'' K in most of the inner crust. The uncertainty in the profile 
of Tc does not affect our results (as can be seen from Fig. |^), since the temperature of the 
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inner crust drops below 10^° K very rapidly (within a timescale much shorter than that for 
the magnetic field evolution in the inner crust). 

3. RESULTS OF NUMERICAL CALCULATIONS AND DISCUSSION 

The three observables are P, P, and P, or equivalently P, Tgd (and -Bsnr/ with 
redundancy), and n. These observables are listed in Table Q together with the parameters 
of our models. Our models aim at fitting these three (four with Bgurf) observed values at 
a given model age tmodei as listed in Table |I]. The age of the Crab pulsar is known, and 
the tmodei is taken to be equal to this age. In the case of PSR 0540-69, the study of the 
associated supernova remnant (SNR 0540-69) indicates an age tsnr ~ 760 yr (Kirshner et 
al. 1989), but the actual age may be somewhat larger or smaller depending on whether the 
remnant's expansion has been accelerated by the pulsar or slowed down by the interaction 
with the interstellar medium. We thus only impose that tmodei ~ ^snr- The age of PSR 
1509-58 is highly uncertain, and tmodei ~ 1,860 yr is an output of our fit for this pulsar. If 
one accepts the association of PSR 1509-58 with the supernova SN 185 (Thorsett 1992), 
then this age is very close to the real age (however, Strom [1994] criticized this association). 

We show in Figure | the temporal evolution of the surface magnetic field Bgurf (strictly 
speaking, its component perpendicular to the spin axis), spin-down luminosity Lg^, spin 
period P, braking index n, and spin-down age Tgd for three different models. The present 
positions (i.e., at age tmodei) on the evolutionary tracks are marked by dots. Figures 5A, 
5B, and 5C show clearly that all observable and derived (such as Bsurf) quantities for the 
Crab, PSR 0540-69, and PSR 1509-58 can be reproduced within our model. Note that 
the calculated values of Bsurf exactly match the corresponding values derived from the 
standard formula for the magnetic dipole losses in which the stellar radius and moment of 
inertia are taken to be equal to those the NS model used for each particular pulsar. 

3.1. Effects of Cooling and Stellar Structure 

Figure |^ shows that with the fast neutrino cooling r^^ exceeds tmodei for a longer time 
than in the case of slow cooling, but the detailed behavior depends also on other quantities 
such as the n and tmodei- The requirement that the model fits the age of the pulsar (besides 
P, Tsd, and n) thus puts strong constraints on the thermal evolution of the NS. In the cases 
of the Crab and PSR 0540-69, such a fit can only be achieved by employing fast neutrino 
cooling, which results in a rapid freezing of the stellar crust and corresponding decrease (see 
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also MP) in the electrical resistivity of the crust. This is clearly seen in Figure after 
t ~ 20 yr, the rate of growth of the surface magnetic field slows down substantially, and the 
magnetic field enters the saturation regime. The Crab pulsar has been claimed several times 
(see, e.g., Tsuruta 1986; Page & Baron 1990) to be a typical example of an NS that follows 
the standard cooling scenario. However, there is actually no real observational evidence in 
favor of these claims (Page 1994, §6.2) and this pulsar may perfectly well have undergone 
fast neutrino cooling as we need here. In the case of PSR 0540-69, the progenitor mass is 
estimated to be around 15-20 M© (Kirshner et al. 1989). The latter implies a high mass 
NS, for which a fast cooling (needed in our analysis to fit the data) is almost inevitable. In 
the case of PSR 1509-58 there is no such age constraint. However, if our model is correct, 
then the fit of P, r^^, and n is evidence in favor of slow neutrino cooling, which, in its turn, 
would argue for a progenitor less massive than that of the Crab [which is estimated to be 
about 9 Mq (see, e.g., Nomoto 1985)]. 

Our conclusion concerning the cooling, fast for the Crab and PSR 0540-69 and slow 
for PSR 1509-58, is quite robust within the framework of our study. However, our analysis 
cannot discriminate between different fast cooling scenarios owing, e.g., to the kaon or pion 
condensates or direct Urea from nucleons and/or hyperons. Note also that our calculations 
are less sensitive to the stellar radius, and the complexity of internal field (the presence of 
higher order multipoles, toroidal component, etc.) may also allow us to relax somewhat the 
constraints on the crust thickness. 



3.2. Predictions for the Crab and SN 1987A 

The detection of thermal emission from the surface of the NS in the Crab can, in 
principle, establish the occurrence of the fast neutrino cooling for this star in the past. The 
new i?05'y4r upper limit (Becker & Ashenbach 1995), based on the absence of detection, is 
already marginally inconsistent with the standard cooling scenarios. Our cooling scenario 
implies an effective temperature well below that in the standard scenario which may 
potentially be discriminated by the future X-ray observatories. Another interesting point 
is that, as one can see in Figure p3, our calculations imply that the spin-down luminosity 
of the Crab has been decreasing only slightly during the past centuries. The comparison 
of the present remnant luminosity with earlier observations may confirm or rule out this 
possibility. We have used an NS model of mass 1.45 Mq for the Crab which is required 
by the need for fast cooling within our TPL EOS. However, this precise value of mass 
should not be taken too seriously, since the critical mass for the kaon condensation is model 
dependent, and the value of a^ms we have adopted in our calculations is arbitrary. Owing 
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to the estimated low mass of the Crab progenitor (Nomoto 1985) one may suggest a lower 
mass for the associated NS. In this case, in terms of NS modeling, wc would simply need to 
increase laam^l to reduce the critical mass for kaon condensation, without affecting seriously 
the maximum mass (which can also be adjusted by changing the underlying baryonic EOS). 
Thus, our estimate of the NS mass for the Crab is not warranted. 

The development of our model was triggered by the announcement of the discovery of 
an optical pulsar in the remnant of SN 1987A. This detection has been neither confirmed 
nor ruled out by Hubble Space Telescope observations (Percival et al. 1995). Our original 
implication ( MP) of the calculations of the field growth for the possible NS in the SN 1987A 
was based on a stellar model with slow neutrino cooling and resulted in a slow increase of 
the pulsar's spin-down luminosity, similar to our present model for PSR 1509-58. This first 
scenario may explain the lack of evidence of powering of the remnant by the pulsar and 
predicts that this powering should be detectable in a few decades or centuries. However, 
our modeling of PSR 0540-58, whose progenitor seems to be a twin of the progenitor of SN 
1987A (Kirshner et al. 1989), argues for fast neutrino cooling in the case of NS 1987A. 
If this very young NS follows similar evolution as that in PSR 0540-58, its spin-down 
luminosity should already be quite high, and it will be growing up by an order of magnitude 
during the next decade or so and may become detectable. If, however, the magnetic field of 
the NS 1987A has been submerged down to deeper layers than in the case of PSR 0540-58, 
then the field growth will be retarded, and the pulsar will be hidden within the remnant for 
a much longer time. The other possibility is, of course, that the NS 1987A did go into a 
black hole (Brown, Bruenn, & Wheeler 1992). This scenario puts strong constraints on the 
maximum NS mass , ~ 1.56 Mq (Bethe & Brown 1995), to which our TPL EOS has been 
adjusted. Also, the scenario of Brown et al. (1992) stipulates that the mass of the NS 
1987A only slightly exceeded this maximum mass. The fact that the NS in PSR 0540-58 
(a probable twin of the NS 1987A) did survive both reinforces their suggestion and may 
imply that it has a mass just below the maximum mass and has undergone a phase of fast 
neutrino cooling. 

4. SUMMARY AND CONCLUSIONS 

We have presented numerical calculations of early magnetic and spin evolutions of 
an NS born with a magnetic fiux initially trapped under the stellar surface. Wc have 
demonstrated explicitly that our calculations may be relevant to the magnetic and spin 
evolutions of the NSs in PSR 0531+21 (Crab), PSR 0540-69, and PSR 1509-58. 

Our principal findings are as follows. 
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1. The braking indices of the Crab pulsar {n = 2.5), PSR 0540-69 {n = 2.0), and PSR 
1509-58 (n = 2.8) can be understood easily in terms of the standard mechanism of magnetic 
dipole braking, provided that the large scale surface magnetic field of the associated NSs 
has been growing up with the rate corresponding to the ohmic diffusion of an internal large 
scale magnetic field, initially submerged under the stellar surface. 

2. We account for the existing difference (~ 300 yr) between the spin-down and real 
ages for the Crab pulsar. Also, according to our calculations, the real age of PSR 0540-69 
is about 890 yr, i.e., 770 yr less than its spin-down age, in agreement with the estimated 
age of the associated SNR 0540-69. Finally, the real age of PSR 1509-58, in contrast to the 
Crab and PSR 0540-69, may exceed its spin-down age by about 300 yr. 

3. In our model the NSs in the Crab, PSR 0540-69, and PSR 1509-58 have initial spin 
periods of, respectively, Pq = 21, 39, and 20 ms, and their initial surface magnetic fields are 
in the range of (1 — 5) x 10^° G. 

4. The spin-down luminosity of the Crab, PSR 0540-69, and PSR 1509-58 reaches its 
maximum at an NS age of ~ 140, ~ 250, and ~ 80 yr, respectively. 

5. Our analysis may sensibly indicate that NSs in the Crab and PSR 0540-69 have thin 
( <200 m) outer crusts and have experienced a phase of fast neutrino cooling. In contrast, 
the NS in PSR 1509-58 probably has a relatively thick (~ 430 m) outer crust, and has 
undergone a slow neutrino cooling. 

To extend our investigation of the early magnetic and spin evolution of an NS to 
other young pulsars, we need to have reliable measurements of their braking indices that 
are not available at the moment. However, the observations recently initiated with the 
ROSAT X-T&j telescope (see, e.g., Aschcnbach, Eggcr, & Triimper 1995) of fine structures 
in young supernova remnants (SNRs) containing pulsars may provide us with more accurate 
estimates of the real ages of young radio pulsars associated with the SNRs. Given such 
estimates, we would be able to proceed farther with a more extensive analysis of the effect 
under consideration in pulsars. 

Finally, a more detailed study should incorporate some other possible mechanisms 
mentioned in the introduction. Spin-up of an NS by early neutrino emission can be 
incorporated in our model since it is a well defined mechanism and requires the same 
ingredients of stellar structure and cooling. The various mechanisms reducing n below 3 
are very model dependent and would increase the number of free parameters. Quadrupolar 
radiation, magnetic and/or gravitational (resulting in n > 3) can reasonably be excluded. 
The crustal "plate" motion is of major uncertainty and is practically impossible to model 
quantitatively. Although plate tectonics is hypothetical, crust cracking ("plate tectonic" 
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activity ?) seems to be taking place in the Crab (Link, Epstein, & Baym 1992; Alpar et al. 
1994), and crustal plates may manifest themselves at the surface of the Geminga NS (Page, 
Shibanov, & Zavlin 1995). 

A.M. thanks the Alexander-von-Humboldt-Stiftung for support. At UNAM, this work 
has been supported by a UNAM-DGAPA grant IN105794 and a Cdtedra Patrimonial from 
Conacyt. The authors also thank an anonymous referee for useful comments, which have 
improved the clarity of the manuscript. 
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Table 1. Characteristics of Pulsars and Parameters of Models 



PSR 


P 


n 


Tsd 


TSNR 


^"rnodel 




Mass 


^sub 


Psub 




(ms) 




(yrs) 


(yrs) 


(yrs) 


(ms ) 


Me 


(m) 


(lO^Og/cm^) 


0540-69 


50.4 


2.01 (2) ^ 


1660 


760 ± 100 


890 


39 


1.50 


85-90 


2.62 3.16 


0531+21 


33.4 


2.509 (1) <^ 


1260 


810 ± 100 


940 


21 


1.45 


95-100 


1.94 - 2.34 


1509-58 


150.2 


2.837 (1) ^ 


1549 


1810 f 


1860 


20 


1.35 


170-175 


1.43 - 1.56 



^Listed are: pulsar spin period P, braking index n, and spin-down age Tsd', associated supernova 
remnant age tsnr', model age tmodel of the NS (as used in our calculations) and initial pulsar spin 
period Pq; theoretical mass of the NS; depth Zgub and density psub of initial submergence of the 
magnetic field below the surface. 

^Nagase et al. 1990; Gouiffes et al. 1992. 
'^Kirshner et al. 1989. 

^^Lyne et al. 1988. 
*^Kaspi et al. 1994. 



^SN 185: Thorsett 1992. 
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Fig. 1. — Mass-radius relationships for the four NS EOSs: FP (Friedman & Pandharipande 
1981), MPA (Miither et al. 1987), PAL (Prakash et al. 1988: see text for the specific 
parameters used), and TPL (Thorsson et al. 1994: see text for the specific parameters 
used). The regions in which M decreases with R are unstable, so that for the TPL EOS 
the branch b-c collapses onto the branch d-c, and the branch a-b is probably not realized 
in nature and also collapses onto c-d. As a result, for the TPL EOS the normal NSs can 
be produced for gravitational masses below 1.37 Mq (point a), while above this mass the 
leap into the regime of formation of nuclear stars (branch c-e with kaon condensate) occurs. 
Point e denotes the maximum mass of 1.575 Mq corresponding to this EOS. 

Fig. 2. — Density profiles for five neutron-nuclear stars from our TPL EOS. At the 
kaon condensation phase transition the density jumps from pi = 1.03 x 10^^ g cm"'^ to 
P2 = 1.76 X 10^^ g cm~'^ and, owing to the enormous amount of kaons present, the proton 
fraction jumps from 9.85% at pi up to 54% at p2, while the electron fraction drops from 
6.14% down to 0.0022%. Despite the large proton fraction and owing to the very low electron 
fraction the direct Urea process is forbidden at densities above p2 and is allowed only when 
p reaches ps = 2.8 x 10^^ g cm~^. See TPL for more details. 

Fig. 3. — Cooling curves for five stellar models with the TPL EOS. The 1.35 Mq star cools 
according to the "standard" scenario, the 1.40 and 1.45 Mq stars undergo fast neutrino 
cooling from the kaon condensate, and the two heaviest stars allow also the direct Urea 
process (see Fig. |^) and cool even faster. The ordinate is the redshifted effective temperature. 



Fig. 4. — Profiles of the internal temperature for two cooling stars with the TPL EOS. The 
profiles for the 1.40, 1.50, and 1.55 Mq stars are similar to the 1.45 Mq case. 

Fig. 5. — Plot of (A) the surface value of the polar magnetic field strength, (B) spin-down 
luminosity, (C) spin period, (D) braking index, and (E) spin-down age vs. the model age of 
the NS. Dots on the curves correspond to the "empirical" and derived data for these pulsars 
(see § § for details). The parameter value Lq in (B) is equal to 2.4 x 10^^, 1.4 x 10^^, and 
4.2 X 10^5 erg s"^ for the Crab, PSR 0540-69, and PSR 1509-58, respectively. The dependence 
Tsd = tmodei IS showu by the dotted line in (E). 
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